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ABSTRACT 

The recent discoveries of luminous quasars at high redshifts imply that black holes more massive than 
a few billion solar masses have been assembled already when the universe was less than a billion years 
old. We show that the existence of these black holes is not surprising in popular hierarchical models 
of structure formation. For example, the black hole needed to power the quasar SDSS 1044-0125 at 
2 = 5.8 can arise naturally from the growth of stellar-mass seeds forming at z > 10, when typical values 
are assumed for the radiative accretion efficiency (~ 0.1), and the bolometric accretion luminosity in 
Eddington units (~ 1). Nevertheless, SDSS 1044-0125 yields a non-trivial constraint on a combination 
of these parameters. Extrapolating our model to future surveys, we derive the highest plausible redshift 
for quasars which are not Icnsed or beamed, as a function of their apparent magnitude. We find that at a 
limiting magnitude of if 20, quasar surveys can yield strong constraints on the growth of supermassive 
black holes out to z ^ 10. 

Subject headings: cosmology: theory - galaxies: formation - quasars: general - black hole physics 



1. INTRODUCTION 

Quasars have long been believed to be powered by the 
accretion of gas onto supermassive black holes (Salpeter 
1964; Zel'dovich 1964; Lynden-BcU 1969). Dynamical 
studies indicate that remnant black holes (BHs) indeed 
reside in the quiescent nuclei of most nearby galaxies (e.g. 
Magorrian et al. 1998; Ferrarese & Merritt 2000; Gebhardt 
et al. 2000, and references therein), and imply that BH 
formation is a generic consequence of galaxy formation. 
Although the activity of bright quasars peaks at z ^ 2.5, 
quasars are known to exist at higher redshifts. Approxi- 
mately 200 quasars have been found so far at z > 4 (the 
majority of which are in the Palomar Sky Survey^), in- 
cluding the record-holding bright quasar SDSS 1044-0125 
at z = 5.80, recently discovered by the Sloan Digital Sky 
Survey (SDSS; Fan et al. 2000). 

Supermassive BHs can grow out of low-mass "seed" 
BHs through accretion or mergers (Rees 1984; Barkana 
& Loeb 2000a). In popular hierarchical models for struc- 
ture formation, the first collapsed gaseous objects have 
low masses. It is natural to identify these galaxies with 
baryonic masses just above the cosmological Jeans mass, 
~ lO"* Mq, as the first sites where seed BHs may form 
(Larson 2000). Here we postulate that supermassive BHs 
are the merger products of individual BHs that grow out of 
these seeds through gas accretion. The natural e-folding 
timescalc for the growth of a single seed can be written as 



where e = Lboi/AfbhC^ is the radiative efficiency for a 
mass accretion rate Mbh, and rj = Lboi/-^Edd is the bolo- 
metric accretion luminosity in Eddington units, Le = 
47rG'AfbhC/ie'Tip/tTT- Here ctt is the Thomson cross-section 
and /ie — 1.15 is the mean atomic weight per electron. 
The growth of a 1O^M0 BH out of a stellar mass seed re- 
quires about ln(lO''A^0/lOMQ) = 18.4 e-folding times or 
~ 7 X 10^(e/0.1)r/^^ yrs. This time is shorter than the age 
of the universe at z = 5.8 if {e/0.1)ri~^ <^ 1. However, the 
time available for the growth of galaxies is shorter than the 
entire age of the universe at the corresponding redshift. 

Turner (1991) recognized that exceptionally bright, 
high-redshift quasars may yield interesting constraints on 
cosmic structure formation. Motivated by the recent dis- 
covery of the SDSS 1044-0125 quasar at z = 5.80 and the 
emergence of a concordance cosmological model for struc- 
ture formation'^, we revisit this problem. In this Letter, 
we derive new constraints from the observed properties of 
known high-redshift quasars, and assess the potential for 
tightening these constraints with future quasar surveys. 

2. COSMOLOGICAL GROWTH OF QUASAR BLACK HOLES 

Our constraints derive from the existence of a BH of a 
particular mass Mbh at redshift z, within a galaxy halo 
of total mass iWhaio- In the following sub-sections we first 
express these masses in terms of observed quantities, and 
then describe the corresponding theoretical modeling of 
the BH growth. 
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■^See tLttp://astro.caltech.edu/''george/z4.qsos for a comprehensive list of known 2 > 4 quasars. 



2.1. Black Hole Mass 



^Throughout this Letter we adopt a ACDM cosmology with Oq 



: 0.3, Oa 
1 



: 0.7, f2b = 0.045, h = 0.7, o-g 



0.9, and n = 1. 
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For the high-redshift quasars of interest, we can only 
infer the BH mass indirectly from the observed quasar lu- 
minosity (although reverberation mapping provides con- 
straints in other cases; see Kaspi et al. 2000). Given 
the luminosity in the observed band, we first apply a 
bolometric correction to obtain the total luminosity Lboi 
(Elvis et al. 1994). For example, the approximate spec- 
tral slope of Fi, oc implies that only ~ 1% of the 
bolometric luminosity of SDSS 1044-0125 is in the SDSS 
z' band. The BH mass then follows from the relation, 
Lbo\ = V^Edd = 1.4 X lO3877(Mbh/M0) ergs-i. In the 
case of SDSS 1044-0125, we infer Mbh = 3.4 x IO^tj-^Mq. 
There are no obvious signs of beaming or lensing in this 
source. However, if present, both of these effects would 
reduce our inferred BH mass or increase our inferred con- 
straint on 7] by correcting the relation between the appar- 
ent and the intrinsic luminosity of the quasar. In §3 we will 
explore the sensitivity of our results to such a correction. 



ponentially by accretion, Afbh(i) = t;xp[At(z)/iacc]-^sccd, 
where At is the time elapsed between the formation time 
of the seed BH and the redshift z. We assume that even- 
tually all massive BHs merge together to form a single 
supermassive BH at the center of the parent halo. As long 
as the BH mergers are completed prior to redshift z, there 
is no need to specify when these mergers took place. The 
mass of the resulting BH in the parent halo is the sum of 
the individual BHs, each of which has grown by a different 
amount, 



Mbh(-2,Mhalo) = Mseed 
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dz' 



■ exp 



At(z,z') 
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where Np^og{z') is the total number of seeded progenitors 
at redshift z' > z, 



2.2. Black Hole Growth 

In order to model the growth of BHs in our adopted 
ACDM cosmology, we rely on the merger history of dark 
matter halos in the extended Press-Schechter (EPS) for- 
malism (Press & Schcchtcr 1974; Bond et al. 1991; Lacey 
& Cole 1993). We first need to determine the mass of 
the halo in which an observed quasar resides. We esti- 
mate Mhaio based on the halo abundance within the EPS 
formalism. The number of halos with mass > Mhaio, 
in a survey of solid angle Af7, probing a redshift range 
Az, is given by 7V(> Mhaio) = AnAz x (dV/dndz) x 
^ dM[dn/dM{z, M)], where dV/dfldz is the cosmo- 
logical volume element, and dri/dM is the number of ha- 
los per halo mass per unit volume. The halo mass for a 
quasar of luminosity L is found by equating the observed 
number of quasars brighter than L to f{z) x iV, where /(z) 
is the fraction of halos hosting active quasars. Because of 
the exponential shape of the EPS mass function for rare 
halos, our inferred halo mass depends only logarithmically 
on f{z), and we conservatively adopt f{z) ~ 1. In the case 
of SDSS 1044-0125, one bright quasar was found within a 
~ 600 deg^ survey area at redshift z = 5.8 ± 0.5; from this 
we infer Mhaio ~ 1-1 x 10"^^ Mq, with a corresponding ve- 
locity dispersion of cr » 470 km s~^. It is interesting that 
the values of Mbh and a we obtain for SDSS 1044-0125 
are in good agreement with their relation measured in lo- 
cal galaxies (Gebhardt et al. 2000; Ferrarese & Merritt 
2000). 

We next compute the growth of the central BH mass 
based on the assembly of its host halo. Given a halo of 
total mass i\/haio at a redshift z, the EPS formalism spec- 
ifies its average merger history at higher redshifts. Every 
branch of this merger tree represents a progenitor of the 
parent halo, whose mass is continuously growing through 
accretion or mergers with other halos. To keep our model 
simple, we assume that every building block of the original 
halo develops a seed BH of mass Mgeed as soon as it ac- 
quires a minimum mass, Mmin, corresponding to a velocity 
dispersion (Tmin- The physical motivation for this choice is 
that radiative cooling and feedback processes which likely 
determine whether or not a massive BH forms in a col- 
lapsed halo, depend directly on a (e.g., Haehnelt et al. 
1998). The mass of each seed BH is assumed to grow ex- 



M r~M- n^'" ,,, dP{z,z',M^,,,,M) 



Here dP{z, z', Mhaio, M) is the number of progenitors in 
the mass range between M and M + dM at a redshift z' 
for a halo whose mass at redshift z is Mhaio [see Lacey 
& Cole 1993, Eq. (2.15)]. Note that we adopt the most 
optimistic assumptions regarding the BH growth; in real- 
ity the growth of each seed BH may be limited by its fuel 
reservoir. Also, the central BHs may not get incorporated 
into a single supermassive BH during galaxy mergers (Hut 
& Rees 1992); however, for exponential growth, the final 
BH mass tends to be dominated anyway by the first seed 
to have formed. 

In summary, our model for the assembly of BHs has 
four free parameters: amin, Mgeed, e, and rj. Although 
the values of these parameters are a'priori uncertain, their 
"natural" choices are as follows: (i) if the seed BHs are 
the remnants of massive stars, then their characteristic 
mass is Mgeed ~ 10 M©; (ii) a necessary requirement for 
star-formation is efficient cooling; in a metal-poor pri- 
mordial gas with no H2 molecules, this implies a mini- 
mum virial temperature of T <; 10*K for progenitor halos 
with seed BHs, corresponding to a minimum velocity dis- 
persion fjmin ^ 10 km s~^; (iii) if the gas accretes to a 
non-rotating BH through a steady thin disk, then the ra- 
diative efficiency, e w 0.06; and (iv) for high fueling rates 
bright quasars would naturally shine close to their limiting 
luminosity, implying 77 « 1. Variations on these parameter 
values will be discussed in § 3 and § 4. 



3. RESULTS 

Given the parameters crmiin Mgeed, and 77, our model 
yields the mass Mbh of the supermassive BH at the center 
of a halo of mass Mhaio at redshift z. In the following, we 
derive constraints on these four parameters, by requiring 
that Mbh equals the value estimated for observed high- 
redshift quasars. 
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No. of Seeded Progenitors Contribution to BH Mass 




Fig. 1. — The assembly history of the black hole (BH) in the 
SDSS 1044-0125 quasar at z = 5.8. The inferred BH mass is 
Mbh = 3.4 X IO^ti-'^Mq and the halo mass is 1.1 X lO" Mq. The 
four panels show, clockwise from top left, the number of seeded pro- 
genitors (those with cr > 50 km s~^); the contribution from different 
redshifts to the final BH mass; the time available for the exponen- 
tial growth of a seed between z and the redshift of 5.8; and the halo 
mass corresponding to 50 km s~^ at each redshift. 



3.1. Illustrative Example 

In order to illustrate the BH growth process in our 
model, we show in Figure |l| the evolution of various quan- 
tities for the SDSS 1044-0125 quasar. In this example, 
we have assumed Mgccd — 10 M0,e — 0.1,77 — 1, and 
o'min = 50 km s^^. As we find by numerical integration 
of equations (Q) and (^), this combination yields the re- 
quired BH mass of Mbh « 3.4 x 10^ M© at z = 5.8. The 
top left panel in Figure m shows the number of progenitors 
of the parent halo with a velocity dispersion > (Tmin- For 
reference, the bottom left panel shows the corresponding 
minimum progenitor mass. As the redshift increases, the 
number of progenitors increases, peaking at z « 9, and 
then decreasing again as the typical progenitors are bro- 
ken up into halos with cr < 50 km s~^. The top right panel 
shows the contribution of progenitors from each redshift to 
the final BH mass, and demonstrates that the bulk of the 
BH mass is contributed by seeds from z « 15. This red- 
shift is considerably higher than the peak at which most 
progenitors form. The increased time available between 
higher redshifts and z = 5.8 (shown explicitly in the bot- 
tom right panel) makes the contribution from the first few 
progenitors dominant. 



3.2. Constraints from SDSS 1044-0125 



Fig. 2. — The minimum value of the typical quasar luminosity in 
units of the Eddington limit {rj = L/L^^^) as a function of the typ- 
ical radiative efficiency of accretion (e = L/rhc^). The constraints 
are derived by requiring that the total BH mass in SDSS 1044- 
0125 will build up to its inferred value of 3.4 X 10^r;~^MQ at 2 = 
5.8. The solid curve shows a fiducial model with (Mgeedi ""min) = 
(10 M0, 10 km s-i); other curves correspond to variations as la- 
beled. 

Next, we search for general combinations of Afgocd, 
77, and (Tmin that yield the required BH mass for the 
SDSS quasar. We solve equation (||) with its left hand 
side set to 3.4 x 1O^77~^M0, and the halo mass set to 
1.1 X 10^^ M0 at z = 5.8, as discussed above. For a given 
pair of (Msocd, o-min), we then find ry as a function of e by 
a Newton-Rhapson method. 

The results of this procedure are displayed in Figure ^in 
a fiducial model with (Mseed, cTmin) = (10 Mq, 10 km s~^) 
(solid line) and its variants. The dotted and short- 
dashed curves show results when (Tmin is either increased 
to 50 km s~^ or decreased to 1 km s~^, respectively. The 
high value applies if the limiting mass for a progenitor 
halo is determined by feedback from the UV background 
following the reionization epoch (Thoul & Weinberg 1996; 
Navarro & Steinmetz 1997; Haiman, Madau & Loeb 1999); 
and the low value applies before reionization if sufficient H2 
exists to allow cooling in small halos (Haiman, Abel & Rees 
2000). The long-dashed curve demonstrates the effect of 
increasing the typical seed mass to 10'^ Mq, correspond- 
ing to remnants of very massive metal-free Population HI 
stars (VMOs; Bond, Arnett & Carr 1984; see also Schnei- 
der et al. 2000). The dot-short-dashed curve describes the 
case where the observed emission from SDSS 1044-0125 is 
beamed into a tenth of its sky, so that its true BH mass 
is reduced by a factor of 10 and the halo abundance is in- 
creased by the same factor. Finally, the dot-long-dashed 
curve shows results for the fiducial model in the case of the 
exceptionally bright z = 4.6 quasar PSS1347-I-4956. We 
have searched through the list of known z > 4 quasars and 
found this object (BH mass of 5 x 10^° Mq, and a host 
halo mass of 2.7 x 10^^ Mq (cr = 570 km s~^, assuming no 
beaming or lensing) to provide the second strongest con- 
straint after SDSS 1044-0125. A handful of other bright 
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z > 4 quasars follow closely behind. 

Figure pleads to two interesting conclusions. First, the 
set of rather standard values, Mscod = 10 M©, (Tmin — 
50 kms"^, e — 0.1, and 77 = 1, yields the required BH 
mass of SDSS 1044-1215. The corresponding progenitor 
halos have cooling times much shorter than the dynami- 
cal time, and could withstand photo-ionization heating by 
the intergalactic UV background (e.g. Thoul & Weinberg 
1996; Navarro & Steinmetz 1997). In such halos there 
is no obvious obstacle to the formation of massive stars, 
which could leave behind the required ~ 10 Mq seed BHs. 
Second, the inferred values of rj are relatively high, sug- 
gesting that BHs must radiate close to their Eddington 
limit. Lower values of 77 are allowed if e ^ 0.05; however, 
such low values of e would not account for the total en- 
ergy output of quasars given the observed mass density 
of their remnants (e.g. Ho & Kormendy 2000; Salucci et 
al. 1999; Fabian 2000). The radiative efficiency may in 
fact approach a value as high as e = 0.42 for a maximally 
rotating Kerr BH. 
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Fig. 3. — The maximum allowed BH mass (lower panel), and 
the corresponding minimum apparent K magnitude (upper panel), 
for a single quasar that could be found within a future survey. 
The solid curves correspond to our set of fiducial parameter val- 
ues, Mgeed = 10 Mq, cTniin = 10 km s~^, € = 0.06, and rj = I; 
the other curves describe variations of this model as labeled. The 
bullets refer to the BH mass and K magnitude inferred for SDSS 
1044-0125. The survey area is assumed to be 10 deg^, except for 
the dot— short— dashed curves which bracket the range between 0.1 
(lower curve on top panel; upper curve in bottom panel) and 10'^ 
(upper curve on top panel; lower curve in bottom panel) deg^. 



3.3. Constraints from Future Quasar Surveys 

In anticipation of future results from the SDSS survey 
as well as from deeper surveys of high-z quasars, we ask; 
what is the highest plausible value of Mbh that should be 
anticipated for a quasar around a redshift z in a given 
area on the sky? To answer this question, we compute the 
mass of the central BH in the most massive dark halo that 
could be found in a representative survey area of 10 deg^ 
and redshift bin Az = 1 around different redshifts. For our 

Small Explorer Mission selected by NASA for study. See a descript 



set of fiducial parameter values we take A/gcod = 10 Mq, 
fmin — 10 km s^^, e = 0.06, and rj = 1. In Figure 
the solid curve in the lower panel shows the maximum 
allowed quasar BH mass in this model as a function of 
redshift. In the upper panel, we show the correspond- 
ing apparent K[AB) magnitude with the corresponding 
K-correction (Elvis et al. 1994). In our fiducial model, 
equation (||) leads to an unphysically large BH mass, in 
excess of the observed halo mass. In Figure || we do not 
allow the BH mass to exceed (57b/f2o)Mhaio- This addi- 
tional constraint appears as the break occurring at z « 9; 
BHs below this redshift are assumed conservatively to have 
consumed all the available gas in their host halo. Figure || 
also shows the maximum allowed BH mass and the corre- 
sponding minimum apparent magnitude, in variants of our 
fiducial model: examples in which (Tmin is either increased 
to 50 km s~^ or decreased to 1 km s~^ (dotted and short- 
dashed curves), and a model in which the typical seed 
mass is increased to 10^ (long-dashed curves). The 
dot-long-dashed curves show the results when we allow 
BHs to consume only 1% of the available gas in its host 
halo. Finally, the dot-short-dashed curves that bracket 
our fiducial model assume survey areas of 0.1 deg^ and 
1000 deg2. 

We find that the survey area enters only logarithmically 
into our constraints. However, the allowed BH mass and 
the resulting apparent magnitude, are strong functions of 
redshift. At its planned sensitivity, the Next Generation 
Space Telescope (NGST) will have a detection threshold 
of ~32 mag in the 1 — 5/im range (with <; 3 hours of 
integration and S/N=10). The assumed survey size of 0.1 
deg^ can be covered by 23 images taken by NGST in a 
total observation time of 2.3 days. The proposed PRIME 
survey** plans to map an area of ~ 10 deg^ to K = 27. 
At a limiting magnitude oi K ^ 2Q (30), future surveys 
can yield strong constraints on the growth of supermassive 
black holes, as they can map out the "brightest quasar" 
envelope shown in Figure ^ out to z 10 (15). 

4. CONCLUSIONS 

Existing data on high-redshift quasars implies that BHs 
as massive as ~ 3 x 10^ Mq were assembled when the 
universe was only a tenth of its present age. Figure 2 
shows that the massive BH inferred for SDSS 1044-0125 
at z = 5.8 can grow in hierarchical galaxy formation mod- 
els with plausible parameter values for the initial seed mass 
(~ 10 Mq), the minimum velocity dispersion of collapsed 
objects that harbor such a seed (^ 50 km s~*), the radia- 
tive efficiency 6%), and the luminosity in Eddington 
units (r^ 1). Figure 3 illustrates the upper envelope of 
plausible luminosity values for high-redshift quasars that 
might be found in future surveys. 

There are several caveats to the constraints we derived. 
First, beaming or lensing may affect the apparent magni- 
tude of some of the brightest and highest redshift quasars 
(Barkana & Loeb 2000b). Second, we have assumed that 
rj and e maintain the same values during the luminous 
quasar phase and the main growth phase of the BH mass. 
It is possible, however, that e is negligibly small during 
the early growth phase of the BH, and high during the 
luminous quasar phase (e.g. Haehnelt & Rees 1993). This 
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would be equivalent to a corresponding increase in the 
seed BH mass in our model. Multiple seeds per halo 
which eventually coalesce are also equivalent to a single 
massive seed with the sum of their masses. Finally, al- 
though a recent cosmological simulation finds good agree- 
ment with the Press-Schechter mass function for the low- 
est halo masses of interest here at z = 10 (Jang-Condell 
& Hernquist 2000); the accuracy of this ansatz still re- 
mains to be tested over a wider range of redshifts and 
halo masses, where there might be systematic deviations 
(Jenkins et al. 2000; Sheth & Tormen 2000). 

Recently, various modifications to the standard ACDM 
model have been proposed because of a potential conflict 
between theory and observations on small spatial scales 
(e.g., Kamionkowski & Liddle 2000; Bode, Ostriker & 
Turok 2000; Barkana, Haiman & Ostriker 2000, and ref- 



erences therein). Our constraints are expected to tighten 
significantly in such models, since they generically reduce 
small scale power and thus eliminate potential hosts of 
seed BHs at high redshifts. Future surveys of high-redshift 
quasars will test these models. 
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